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1 Introduction 

Non-leptonic B decays have great affluence of phenomenological applications. Especially, 
they provide a plenty of fascinating CP asymmetries in the standard model (SM), where we 
have a number of chances to test the SM by searching discrepancies between the theoretical 
predictions and experimental measurements. The measurements of CP asymmetries have 
been performed well in various B factories such as Belle and BaBar and also in Tevatron. 
Furthermore, upcoming B factories such as LHC-b, Super-B and Belle-2 will shed more 
light on CP asymmetries by accumulating much larger amount of data. 

For a long time, B — )■ Kn decays have been a continuously engrossing issue due to 
that some puzzling behavior occurs in their CP asymmetry measurements: The difference 
between Acp{Tr^K ) and Acp{tt^K ) measured by Belle, BaBar and CDF collaborations 
is unexpectedly larger than the theoretical estimates [1-4] in the SM. Recent brand-new re- 
sult from the LHC-b experiment also supports the behavior [5]. There have been numerous 
new physics analyses on the issue model-independently or model-dependently However, 
before we explore new physics effects beyond the SM to resolve the puzzle, it is strongly 
required that one should perform more precise higher order corrections to the observables 
within the SM. 

The QCD factorization (QCDF) is a theoretical framework for systematically calcu- 
lating hadronic matrix element of weak decays of B meson [6]. It has been extensively 
applied to non-leptonic B decays [7-9]. The QCDF formalism for hadronic matrix element 



of non-leptonic B decays is described by 

{M^M2\Qi\B) = Y,Pf^'''"i^l)fM2 j duTl^{u)^M,{u) + (Mi ^ Ah) 

+ fBfMjM, I didudvTl\i,U,v)^B{i)^MM^MM , (LI) 

where F ^ ^ is B ^ Mi form-factor, /a/- is decay constant of meson Mj, and ^Mi{u) is 
light-cone distribution ampHtude of meson Mi with parton momentum fraction u. The for- 
mahsm is expressed by convolution of hard-scattering kernel T ' with meson distribution 
amplitude at the leading power of h.QCD/fnh- The hard-scattering kernel is separated into 
the hard-scattering form- factor term T^ and the hard-scattering spectator term T^^ . The 
very first paper of QCDF calculated next-to-leading order (NLO) correction to i? — t- vrvr 
decays [6]. The imaginary part of the decay amplitude, which causes strong CP phase, 
arises first at the order Ug from hard-scattering contribution between decaying quarks of 
different mesons. Therefore, the next-to-next-to leading order (NNLO) correction is of 
great importance in perturbative expansion in order to provide reliable prediction for CP 
asymmetries. 

The NNLO correction to tree amplitude for hard-scattering form-factor term has been 
calculated first for imaginary part [10] and later for real part [11], where the order a1 
contribution is quite significant compared to the order Us contribution. Especially, the 
color-suppressed amplitude is very sensitive to the order a1 contribution. The result is 
also confirmed by ref. [12]. As for the hard-scattering spectator term, the one-loop order 
Og correction has been calculated for tree amplitude [13-15] and for penguin amplitude 
[16, 17]. However, the order a^ penguin correction for the hard-scattering form- factor 
term has not yet been calculated. This correction is more important for B — t- Ktt decays 
because the tree amplitudes are CKM-suppressed and the decays are penguin-dominant. 
Especially, the order a^ penguin correction is highly required for theoretical estimate of 
CP asymmetries since it is the first correction in perturbative expansion. 

Motivated by current status of higher order correction for B to light meson decays, here 
we compute the order a^ one-loop contribution (NLO) for the magnetic penguin operator 
as the first step toward the complete order a^ correction of penguin amplitude. In section 
2, we describe the formalism for the calculation. In order to regularize ultra-violet (UV) 
divergences and infra-red (IR) divergences in the NLO calculation, we use dimensional 
regularization where the space-time dimension d is analytically continued to d = 4 — 2e. 
Especially, we choose naive dimensional regularization (NDR) scheme for the prescription 
of 75. We use 't Hooft-Feynman gauge (^ = 1) throughout the calculation. In section 3, we 
explicitly show that the soft and collinear divergences are canceled out in the calculation. 
The result of calculation is expressed in section 4 in complete analytic forms. We relegate 
our analytic results of the master integrals to appendix. We apply our result to -B — > Kit 
decays and provide numerical values of the NLO correction for magnetic penguin operator 
compared with other contributions. In section 5, we summarize and conclude. 



uq 



q 

^- 



-uq 



Figure 1. Order a^ (LO) diagram for magnetic penguin operator with momentum configuration. 
M = 1 — u is understood. The square box denotes Qgg insertion. 

2 Formalism 



The effective amplitude of non-leptonic weak decays of B meson, such as i? — ^ M1M2, is 
expressed [18] by 

^-ff = ^E^'C',(^)(MiAf2|Q.|i3)(/i) (2.1) 

with effective operators Qi, i = 1, 2, ..., 10, 77, 8g. Gp is Fermi constant and A* is a factor 
for CKM matrix elements corresponding to operator Qi. The magnetic penguin operator 
Qgg) on which we are focusing in this work, is defined by 



Qsg = -^^mhsa^y{l + 75)G^''6, 



(2.2) 



where C^*^ is gluonic field strength tensor contracted by SU(3) generators. We do not 
consider the magnetic 7-penguin operator Q^^ because the contribution is much suppressed 
by the fine-structure constant. The Wilson coefficients Ci{^) are to be calculated near weak 
boson mass scale M^y where we can safely avoid large logarithms in perturbative expansion. 
Then, Ci{^) are evolved into ttij, scale through renormalization group equation associated 
with anomalous dimension matrix. In principle, the scale dependence of C.j(/.f) should be 
canceled by the scale dependence of the full calculation of matrix element {MiM2\Qi\B){^) 
up to the order of truncation of perturbative expansion. 

The calculation of hadronic matrix element can be handled within the QCDF frame- 
work. The formalism is expressed in eq. (1.1). In this framework soft gluon exchange 
between decaying quarks is power suppressed in heavy quark limit. The dominant hard 
gluon scattering contribution is absorbed into hard-scattering kernels T ■ . The matrix 
element is described by convolution of hard-scattering kernel with meson distribution am- 
plitude. Basically, T starts at order unity in perturbative expansion, while T starts 
at order Ug- It should be noted that the contribution of magnetic penguin operator in 
T starts at order as- In order word, leading order (LO) diagram for magnetic penguin 
operator is order Us as shown in figure 1. This property is same with all other penguin 
amplitudes. In the figure 1, spectator quark line is suppressed. Left horizontal line denotes 
incoming b quark with momentum p and right horizontal line is outgoing quark with mo- 
mentum p' which forms a bound state of a meson with spectator anti-quark. The upper 
side lines represent the quark and anti-quark that make a bound state of another meson 



with momentum q. u and u (equal to 1 — u) is parton momentum fraction of each quark 
in the upper side meson. The kinematic relations among the momentum four vectors are 
as follow: 







P^ = ml, q^ = 0, 


p" = o, e = 


uml, 


■1-- 


= p-p' -- 


/ / ml 

= p ■q = q-q = ^, 


y-.' = ";l 


p-q' 



ii±|H . (2.3) 

In order to obtain hard-scattering kernel T (superscript / is suppressed) for magnetic 
penguin operator, we use the QCDF formula which is schematically described by 

(Q8s)rcn = i^-r®$M, (2.4) 

where the symbol ® denotes convolution with meson distribution amplitude. {Q8g)rcn is 
renormalized matrix eleinent which is related with its bare quantity by 

{Q8g)vcn = ZmZgZqZ^ ^88(Q8g)0) (2-5) 

where Zgs is operator renormalization constant for Qsg, and Zm, Zg, Zg, Zq are the 
mass, coupling, quark field and gluon field renormalization constant, respectively. The 
one- loop calculation of Zgg is shown in ref. [19]. We summarize each value of one- loop 
renormalization constant: 



Zm = -3Cf , 
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nf, 



^p) = -c^. 



<? 



Z'i/^ = 8Cf - 2N, , (2.6) 

where the following power expansion is considered 

oo j .. 

^. = i + EE(S)'?#"- (") 

i=i fc=i 

Each quantity in the factorization formula eq. (2.4) can be power expanded. Since {Qsg)ren 
as well as T start at order a^, we describe 

(Q%>,.« = (^) {{Q>,& + (^) (<?%>!;!, + 0(af)) , 

From the power expansion of factorization formula eq. (2.4), we obtain 

(Q8,)(°) = r(0) , (2.9) 

(Q8^)(i) = FWr(i) ® ci>W + F(i)T(0) ® ci>W + F(0)r(0) » ci>W . (2.10) 
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Figure 2. NLO diagrams for magnetic penguin operator. Tlie momentum configurations for all 
the diagrams are same as in figure 1. The square box denotes Qgg insertion. 



Those eqs. (2.9,2.10) are the master formula for the calculating hard-scattering kernels. In 
those formula, all the quantities are renormalized while they might have IR divergences: 
All the IR divergences should be canceled in those equations. In the next section it will be 
explicitly shown that the soft and collinear divergences are canceled out in those equations. 
The NLO diagrams for magnetic penguin operator are shown in figure 2. 

In order to compensate scale dependence of NLO hadronic matrix element, we have 
to employ next-to-leading logarithmic (NLL) Wilson coefficients for the Qsg- We refer to 
refs. [20, 21] for NLL effective Wilson coefficient Cg^ at low energy scale. The light-cone 
projection operator of a light pseudoscalar meson in momentum space is applied to upper 



side quark and anti-quark spinors. The projection operator in leading twist is given by 



M, 



a/3 



4iV, 



IsI^^^mN, 



(2.11) 



where a, ^ are spinor indices of upper side quark and anti-quark spinors. The distribu- 
tion amphtude of meson M denoted by $m [u) is conventionally expanded in Gegenbauer 
polynomials: 



^]\l{u,i^i) = 6u{l -u) 



l + ^at'{^)C^^/'\2u-l) 



n=\ 



(2.12) 



.(3/2) 



a^ (fj.) are Gegenbauer moments and Cn are Gegenbauer polynomials of order |. We 
truncate the expansion by n = 2. In order to consistently compensate NLO hadronic 
matrix element, we also use NLL evolution equation for Gegenbauer moments which are 
obtained in [22-24]: 
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where the anomalous dimension matrix elements are given by 7J 
one-loop order, and 7^ = -^, 4'^ = '2°°° 
each number. 
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(2.13) 
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2^ for two-loop order. We set nf = 5 for 



3 Soft, Collinear and UV divergence 



3.1 Soft and collinear divergence cancelation 

The soft divergences can occur when the loop momentum / goes to 0. The collinear diver- 
gences can arise when I becomes collinear to the momentum of the massless quark q or p'. 
All these soft and collinear divergences should be canceled out in order to verify that the 
QCDF framework is valid. In this subsection we explicitly show that the soft and collinear 
divergences of the order a^ correction for magnetic penguin operator are canceled out in 
the QCDF framework. 

First, we examine the soft divergence cancelation. We set the loop momentum / ~ A 
where A is approximately which can be used for investigating the degree of soft divergence. 
In this way we can easily check the soft divergence by counting the power of A in the integral. 
We note that d'^l ~ A^. If the power of A is smaller than or equal to 0, the integration has 
soft divergence. 

It turns out that the diagrams (le), (If), (2) and (5a), (5b), •••, (7d) have no soft 
divergences. The soft divergence of the diagram (3) vanishes due to the on-shell condition. 



The discussion about the diagram (4) will be done in the later part of this section. We 
now consider diagrams (la) and (lb). Except overall factor including color factor (if not 
mentioned, we imply this suppression below), their integrand is 

(la) + (16) 






{I + q')^l + uq)^ 



(3.1) 



where we suppress the spinor indices. We use square brackets implying separate spinor 
indices. T^, is a Dirac structure coming from Qsg insertion which is defined by 

r^, = <T'-g:,(l + 75). (3.2) 

After I goes to and doing some Dirac algebra, we easily find the terms in the round 
brackets cancel each other. Similarly, one can find that the soft divergences from diagrams 
(Ic) and (Id) are also canceled each other. Then, we reach that all the soft divergences are 
safely canceled. 

For the collinear divergence, we first consider the case where loop momentum I becomes 
collinear to the momentum q of upper side quark lines. We typically decompose I as follows, 

l = aq + f3q + ll, (3.3) 

where q = {E,0,0,E), q = {E,0,0,—E) and l± = {0,li,l2,0)- E is energy of upper side 
decaying meson which is approximately E ~ mt/2. Because / is collinear to q we set the 
order of each parameter such that a ~ 1, /3 ~ A^ and l'^ ~ A^m^, where A ~ AqcD/w-fe- We 
note that d'^l = l/2dlqdlqdl\ ~ A^m^ and P ~ A^m|. Then we count the power of A in the 
integral which implies the degree of collinear divergence. Similar to the soft divergence, if 
the power of A is less then or equal to 0, the integral has collinear divergence. One may 
simply check the collinear divergences from the fact that the gluon, which is connected any 
one of the external massless quark line, potentially has collinear divergence. 

It is simple to find that the diagrams (2), (5a), (5b), •••, (7c) have no collinear di- 
vergence. It is worth showing the integrand of diagram (7d) for studying the collinear 
divergence: 

(„) = --l!%^, (3.4) 

where V^^ is the Dirac structure from the qqgg vertex of Qgg which is defined by 

f'^- = c7'^-(l + 75). (3.5) 

After substituting I for aq in the numerator and using some Dirac algebra and on-shell 
condition, we find that the numerator vanishes. Now we consider again diagram (la) 
including color factor: 



Cf\ 1 



<'"' = ( - iB^HnfT^^^-^"^^ + "^'^«'7T73^^'"ih'-: 



(3.6) 



After we substitute I for aq and using ijiad +f — rrih) = 1/a from the on-shell condition, 
we arrive at the cohinear divergence term ^(la) '■ 

where P{u) is LO contribution of magnetic penguin operator defined by P{u) = Cp[T'^,]['y^]/{um'^). 
Similarly, we can find the collinear divergences of diagrams (Id) and (le): 

■V) = (2S; + '^-)p(7TT^(^)<-^<" + °''' P-^> 

From the definition of P{u)^ one can easily find the relation 

P(u + a) = - -P('u), (3.10) 

[u — a) 

where we used on-shell condition in order to get [F^,] = \^^i]- Using this equation we can 
find that the collinear divergences proportional to the color factor l/{2Nc) in X/;^^),^^^) 
and X(ie) are exactly canceled each other. The remaining collinear divergence term is 
proportional to Cp-, and represented by 

^(ia)+(i.)+(i.) = ^^^2(7^ (^) (^W - nu + a)) . (3.11) 

We also consider the collinear divergences of diagrams (lb), (Ic) and (If), and the calcu- 
lation is similar to previous one. The result is 



V / CK + li \ 



'a + w 

Kl6)+(lc)+(l/) = O^p^T^^l^^ 



For the diagram (3), if we consider the projection operator eq. (2.11), we get the following 
collinear divergence term: 



21] 
P{1 + uqY{l - uqy 



^(3) = CFTTTTT-Zyhl ^Zv^^i^ + ") • (3-13) 



After we change a into —a in eq. (3.12) and take apart the propagators in eq. (3.13), we 
obtain total collinear divergence term combining eqs. (3.11), (3.12) and (3.13): 

^ /2(a + n) 1 2{u-a) 1 \,„, , „, ,, 

As this equation is same to eq. (193) in ref. [7] except T(n) is changed into P{u), this 
collinear divergence is to be canceled by one-loop correction of meson distribution ampli- 
tude. To be self-contained, here we show again how this cancelation take place. There is 



missing coUinear divergence that conies from the self energy diagrams of upper side quark 
fields. We find that this additional collinear divergence is given by 

X.a^-C.(<^,,^-i^,,i^)p(„). (3..) 

We add Xadd to Xtot and integrate them over (3 and l± in order to express them as a 
convolution in a. Using Cauchy's theorem for integration over /3 and 

^ = 21n^ 

we can find that Xtot + -'^add is expressed by 

^tot + Xadd = Cf— In ^^ / dwP{w)V{w, u) . (3.17) 

TT ^J'm Jo 

V{w,u) is the Efremov-Radyushkin-Brodsky-Lepage (ERBL) kernel [25, 26] defined by 



^ _21n^, (3.16) 



V{w,u) 



e{u -w)-(l + -^— ] +6{w-u)^{l 



u \ u — w J u \ u — w 

with the following definition 



(3.18) 



[f{w, u)]+ = f{w, u) - 5{w -u) I dvf{v, u) . (3.19) 

Jo 

We note that the NLO contribution of light meson distribution amplitude is expressed by 

$W(w;) = Cf— In ^^ / duV{w, u)<^fhu) . (3.20) 

vr fiiR Jo 

Therefore, after the convolution with meson distribution amplitude, the total collinear 
divergence Xtot + -'^add in the matrix element is exactly canceled by NLO correction to the 
meson distribution amplitude in eq. (2.10). 

Finally, we consider the collinear divergence that arises when the loop momentum / 
is collinear to q' . Diagrams (Ic), (Id), (If), (4) and (7d) have this collinear divergence. 
Without showing the details, it is found that the collinear divergences of diagrams (Ic) and 
(Id) are canceled each other and so do the diagrams (If) and (7d). The diagram (4) has 
both soft and collinear divergences. It is straightforward to show that the soft and collinear 
divergences of diagram (4) are exactly canceled by second term of eq. (2.10) which is 
contributed by NLO form-factor correction diagram shown in figure 3. The UV divergence 
of NLO from-factor correction should be canceled by form-factor renormalization. 

3.2 Renormalization - UV divergence cancelation 

Now we add gluon self-energy diagrams as well as its counter term as shown in figure 4. 
The third diagram is the counter term for UV cancelation of diagrams (Ic) and (If). Other 
diagrams that cause UV divergence are diagrams (la), (le), (2), (5a), (5b), (6a-6d) and (7a- 
7c). These diagrams are same with those for calculating operator renormalization constant 
Zgg except right side quark current. Therefore, all these UV divergences are naturally 
canceled by operator renormalization equation (2.5). 



Figure 3. NLO form- factor correction diagram. 
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Figure 4. Gluon self-energy diagrams and the counter terms. 



4 Result and application 



For the calculation, we first reduce the tensor integrals into scalar Feynman integrals. In 
order to reduce scalar Feynman integrals into the master integrals, we use the Laporta's 
algorithm [27] with in-house Mathematica code which facilitates several reduction methods 
such as Passarino-Veltman reduction [28], integration-by-part method [29, 30] and Lorentz- 
invariance method [31]. We mainly use the Mellin-Barnes representation to compute each 
master integral. Analytic results of the master integrals are shown in appendix. We 
explicitly separate the 1/e and 1/e^ terms into IR divergence part and UV divergence part 
by investigating the origin of divergences for each scalar Feynman integral. We confirm 
that all the IR divergences and UV divergences are canceled separately. Here we show 
the result of Tq and Tq (we attach subscription Qg in order to denote contribution of 



magnetic penguin operator) except the color factor Cp/Nc- The LO result Tq 
same as in refs. [6, 8]. 



-2/u is 
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The NLO result of Ti reads as 

.„(c.((.„--)M.)-f-i.2) 

+ ^((2u-^) 1.1 (S) - i|^ - - - 2u 1.1(1.) - 1 

^ f f 2u \ , ■? , ^ /8u 4 \,,^ 88n 

+^4(^-2-)ln^(^)+(y + - + 2)ln(.)- — +2 

+^ ( f- + l-3nVi('")- f2 + -V2(^^)-6nLi2(u)-2nLi2(u) 
Nc \\ u J \ uJ 

+ ( — 2 + u J In (u) — n In (u) — Au In('u) ln(u) 

+ (— + --l)ln(^) + 2nln(n) + ^ ^J 

+ ("-/ - 2)^8(0 -ie,u) +gs{sc- ie,u) +58(1,^), (4.1) 

where Sc = ml/ ml . The definition oi g%{s,u) is 

.8(.,n)= — + -n J,(-)+^--.ln(.) + — . (4.2) 



The function Ji {u) and J2 (ti) are defined by 

l+jy(u) 



Li2(n) 7r2 /■! InO-^^O^O) ....^ 

u Qu Jq (1-0(1-^0 



where 



\/4 - u - \/^ 
y(n) = ^==—^= . (4.4) 



\/4 — U + -y/^U 

In the case s — t- 0, we find 

20u 4 4 

5r8(0-ie,'u) = — -■uln(u) + i7r-M. (4.5) 

Even though J2{u) has integral form due to a non-trivial master integral, after convolution 
with meson distribution amplitude we could obtain complete analytic result of the hard- 
scattering amplitude. It should be emphasized that uu Tq has no singularity except log- 
singularity as w — ?• nor u — )• 1, so that there is no end-point singularity for the convolution 
with meson distribution amplitude at leading twist. 
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After convoluting the hard scattering kernels with meson distribution amplitude at 
leading twist, we are left with following formula: 

Ps% = /' du 4>A/(n) = -6(1 + af + af ) , (4.6) 

Psii ^ I' duT^yMiu) 

.(f6.C.(|)-..Cl.(|)-.Cl3(|).^.^ 

- 49V3vr - ^) af + (-96vrCl,(|) + 192Cl3(f ^ 

- '-^^ - m.' + 257^/3. + H^) af + (n, - 2)0^(0 - ^e) 

6 io / 

^G-is. - i.) + Gf (1) + (-^ - ^ - 36) in t ^\ 



3 3/ Vm;,/ 

+,. ((52.= - 1^) „;'= + (5f? - 212.=) „,"'» - ^ + 9) . (4.7, 

The function Gg^(s) is defined by 

G8^(s) = -24^2 ln2 {y{s)) + (4 - 40s) ([^||jy) 1^ (^s)) - 104s - 4 ln(s) + y 



-af ( 24(8s - 9)s' In^ (y(s)) - 4(48s' + 34s - 1) ( ^ _ V,' ] In (y(s)) 



- 192s^ - 440s - 41n(s) + 18 
-af ( - 48 (45s2 - 40s + 18) s^ In^ (y(s)) 
+ 4 (540s3 - 390s2 - 70s + l) f ^^^44^ In (y(«)) 

Vi-y(^)/ 

+ 2160s^ - 1740s2 - 1064s - 41n(s) + 21 ) , (4. 



where 



m = ^^"' . (4.9) 

In order to obtain the above equation, we use inverse binomial summation formula [32]. 
The Clausen function Cln(0) is defined in terms of poly-logarithm function [33] 

Ci (r)-i^ (^^" (^"'') " ^^" ('''')) " ^' ''''™ (4 10) 

In eq. (4.7), we explicitly show the imaginary part in the last line. The other sources of 
imaginary part are Gg'^(0 — ie) and Gg^(sc — ie). These imaginary parts arise first at NLO 
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(order a^) for the magnetic penguin operator which causes strong CP phase in non-leptonic 
B decays. 

Now we apply our result to B ^ Kit decays. They are well studied in ref. [8] within 
the QCDF framework up to NLO. Here we follow the notation and the formula in ref. [8] 
for comparing our order a^ result with their NLO one. The B — )■ Kit decay amplitudes 
A{B —7- ttK) without annihilation amplitudes are expressed by 



A{B' -^ TT-lC) = Xp 



P P \ I Kip P 



A, 



nK, 



3 

Ka2 + Ap-(-a7 + ag) 



-y/2A{B- -^ TT^K^) = [A„ai + Ap« + a%) + Xp r^{al + a^)] A^k 

Aktt, 

-A(^ ^ TT+K-) = [Xuai + Xp{al + <„) + Xp r^ {a^ + a^)] A^k , 
V2A(B° -^ TT^lf) = AiB~ -^ TT^K°) + V2A{B^ -^ tt'^K") 

-A(B° ^TT+K~). (4.11) 

Here, Xp represents CKM factor defined by Xp = VpbVpg, and the summation with p = u,c 
is implicitly considered. It should be emphasized that the terms with A„ is highly CKM- 
suppressed as estimated by |VufeV„*s/V"cfel4s| ~ 0.02. At^k and Axn contain all the hadronic 
parameters such as form- factor, meson decay constant and mass factor. The dimensionless 
coefficients Oj represent the hard-scattering contribution of each operator combined with 
its Wilson coefficient, r^ is chiral-enhancement factor defined by 

-x = -^-\-^ ' (4-12) 

where q = u for charged kaon and q = d for neutral kaon. It is known that a^ and r^Og make 
dominant contribution in i? — t- Ktt decays among the aj's. Both come from the contribution 
of penguin operators. The magnetic penguin operator contributes to both a^ and r^a'^, 
but with different twist order of light-cone projection operator of pseudoscalar meson: 
leading twist for a1 and twist-3 for r^a^. Since, up to order as, the twist-3 contribution of 
magnetic penguin operator is quite suppressed than the leading twist contribution, here we 
only consider NLO contribution of magnetic penguin operator to a^ which requires leading 
twist projection. 

Without considering hard spectator scattering term (the term specified o^jj), we sep- 
arate coefficient a^ j (subscript I denotes the hard-scattering form- factor term) into three 
terms in order to compare the contributions of different sources: 

c c,V , c,P I c,Pg r A -io\ 

04,1 = 04,1 + 04'^ + 04^ , (4.13) 

where a^' j denotes vertex correction term, a^' j represents penguin correction term and a^' j ^ 
is reserved for magnetic penguin operator contribution. We note that the three quantities 
are separately gauge- invariant, since they come from different gauge-invariant effective 
operators. Some care is needed for renormalization scheme dependence in each term. We 
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Table 1. Input parameters. 



Parameters 
Values 



Ag^ m,{mt) m,{mt) Qf(lGeV) af(lGeV) 

225 MeV 4.2 GeV 1.3 ± 0.2 GeV 0.3 ±0.3 0.1 ±0.3 




see that the vertex correction starts from order unity while the penguin correction as well 
as the magnetic penguin correction start at Og, the vertex correction term a^'j makes 
dominant contribution. 

We show again the NLO expressions of each coefficient from ref. [8] including our result 
of NLO magnetic penguin operator contribution 

<J = C4 + ^ (l + Cf^ (l21n !^ - 18 + ^' dug{u)^K{n)^'^ , 



+(C4 + C6)f ^In^ - {uf - 2)Gk{0) - Gk{sc) - Gk{1 



c,Ps _ ^esCpOs^ /p(0) (a^\ pil) \ ,. .. 

«4,I - ^8, ^^ 4^ {^S,K + {^J \k) ' (4.14) 

where the loop-effect functions g{u), Gk{s) are defined in ref. [8], and the P^ ^ can be 
read off from eqs. (4.6) and (4.7). 

For the numerical analysis, we compare the NLO contribution of magnetic penguin 
operator with its LO contribution and also with other penguin contribution and vertex 
contribution up to the NLO. We summarize input parameter values in table 1. We use 
two-loop running coupling constant and running quark mass with provided parameters. 
For the Gegenbauer moments, we use the values at the fixed scale with conservatively 
chosen errors that are consistent with several QCD sum rule results [34, 35]. With this 
preparation, we compute the values for a^'j ^ at the LO and NLO with three different scale 
values. The result is described in table 2. For comparison, we also show numerical values 
of the a^' J and a^' j at NLO. The errors propagated from the uncertainties of Gegenbauer 
moments and charm quark mass are given in the round and square brackets respectively. It 
should be noted that the contribution of magnetic penguin operator has opposite sign from 
the vertex and penguin correction, and compensates the other contribution. Specifically, 
the real value of LO magnetic penguin contribution 04 j* (order as) is comparable to that 
of NLO penguin contribution a^' j (order a^) and both cancel each other at // = rrih. 
As can be seen, the NLO contribution of magnetic penguin operator is significant and can 
considerably reduce the absolute value of total coefficient a^ j. Especially, for the imaginary 
part, a4j*(NL0) contribution is so strong that the contribution significantly reduces the 
total value. We emphasize that the large NLO contribution is mainly from the order a^ 
QCD effect of the hard scattering term. We also comment that the contribution from the 
order q^ hard spectator scattering with QCD penguin contraction is very small, as can be 
seen in ref. [16]. 
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Table 2. Numerical result of coefficients a^'j" up to LO and NLO compared with NLO values 
of 04 J and a^'j at three different scale values. The numbers in round brackets indicate maximal 
variation due to the uncertainty of Gegenbauer moments while the numbers in square brackets 
represent the variation from the charm mass uncertainty. 







Real part 






Imaginary part 




M 


mb/2 


rub 


2mb 


mb/2 


m-b 


2mb 


"■ij 


-0.046(1) 


-0.033 


-0.024 


-0.003(1) 


-0.001(1) 


-0.001 


C,P 

«4,I 


-0.003(2)[1] 


-0.009(1)[1] 


-0.014(1)[1] 


-0.001(3) [4] 


-0.003(2) [3] 


-0.003(2) [3] 


al'p (LO) 


0.014(2) 


0.009(2) 


0.007(1) 











<;f« (NLO) 


0.019(4) 


0.014(3) 


0.010(2) 


0.007(2) 


0.004(1) 


0.002(1) 



5 Conclusions 

We computed the one-loop (order a^) contribution of magnetic penguin operator for B 
decays to light meson within the QCDF framework. It is explicitly shown that all the soft 
and collinear divergences are canceled out. All the results are expressed in complete analytic 
forms so that they can be easily applied to various non-leptonic B decays. It turns out 
that the NLO contribution of magnetic penguin operator is quite significant especially for 
imaginary part in B — )• Kit decays. Therefore, the missing part of the order a^ correction 
to the hard-scattering form-factor term for penguin amplitude is strongly required for 
precise estimates of CP asymmetries. It is interesting that the order a^ contribution of 
magnetic penguin operator considerably reduces the absolute value of penguin amplitude 
coefficient 04. Combining future analysis for the order a^ penguin correction together with 
the current result of order a^ vertex correction and hard-spectator scattering corrections, 
we could reach at the complete order a^ accuracy for B — >■ Kir decays at the leading power 
of AQCD/rrib. 



A Appendix: Analytic formulae of the master integrals 



In this section, we summarize the result on the computation of the master integrals. We 
end up with 8 master integrals which are displayed in figure 5. The calculation is based 
on MS renormalization scheme. Our convention for the master integral with propagators 
Vi,V2,.:,Vnis 



(fie 



7e/2\4- 



d'^l 



1 



ITT' 



d/2 p^p^...p^ 



(A.l) 
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uq 1 \ uq 
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/ 

-P 


"'t^~-r''' 






1 uq 






hi 


hi 



uq 



Figure 5. Master Integrals. Solid line represents 'massive' for propagator whose mass is explicitly 
shown in the diagram, and non-zero virtuality for external line. Dashed line implies massless 
propagator or zero virtuality. All the momenta are understood to be ingoing. 



h 
hi 



--'(^)'(«^nr(-i + e; 



2 



+ iO (eT^^) 



r(i-eme) 



h2{q'<m')= (^ 



m" 



r(2 - 2e) 
x + l)\n{x + l) 1 



+ [- + ! \\n^x+l) 



/22(r =0) 
hz 



X 

where x 

2 



+ - + e^- + l)Li,(-x) 
X 12 






m? 



^ 



m^ 



^ + l + e (l + ^)+0(e' 



]- + 2 + \±y^Hy{-x)) + 0{e) 
e l-y{-x) 



hi — n I n 



^,'^e 



2 \ 1 



l-(^{ln{u)-znf + 7,']+0{e) 



h"? — 7T 



1 //.^ 



2 \ 2 



Li2(u) vr 



u 



Jn(l-t2^(l-0) 



6u i ^^-(l-e)(l-nO 



+ 0(e) 



/di = -^ 



ml 



^ + iO) — 
mf^ J uu 



e^ 



2 In iuu) 



£ 



-2Li2(n) - 2Li2 (n) + In^ {-) +^ + 0{e) 

\uJ o 



-(42 — T I n I — 



4 \ 2 , 



3 In(n) — ivr 13 vr^ ,^, , 



(A.2) 
(A.3) 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 



(A.9) 
(A.IO) 
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We note that maintaining integral form in I32 is convenient for convolution with meson 
distribution amplitude. 
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